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A mid-infrared spectroscopic study of submillimeter galaxies: 
luminous starbursts at high redshift 

E. Valiante 1 , D. Lutz 1 , E. Sturm 1 , R. Genzel 1 , L.J. Tacconi 1 , M.D. Lehnert 1 , A.J. Baker 2 

ABSTRACT 

We present rest frame mid- infrared spectroscopy of a sample of 13 submillime- 
ter galaxies, obtained using the Infrared Spectrograph (IRS) on board the Spitzer 
Space Telescope. The sample includes exclusively bright objects from blank fields 
and cluster lens assisted surveys that have accurate interferometric positions. We 
find that the majority of spectra are well fitted by a starburst template or by 
the superposition of PAH emission features and a weak mid-infrared continuum, 
the latter a tracer of Active Galactic Nuclei (including Compton-thick ones). We 
obtain mid-infrared spectroscopic redshifts for all nine sources detected with IRS. 
For three of them the redshifts were previously unknown. The median value of 
the redshift distribution is z ~ 2.8 if we assume that the four IRS non-detections 
are at high redshift. The median for the IRS detections alone is z ~ 2.7. Placing 
the IRS non-detections at similar redshift would require rest frame mid-IR obscu- 
ration larger than is seen in local ULIRGs. The rest frame mid-infrared spectra 
and mid- to far-infrared spectral energy distributions are consistent with those 
of local ultraluminous infrared galaxies, but scaled-up further in luminosity. The 
mid-infrared spectra support the scenario that submillimeter galaxies are sites 
of extreme star formation, rather than X-ray-obscured AGN, and represent a 
critical phase in the formation of massive galaxies. 

Subject headings: galaxies: starburst, galaxies: active, galaxies: distances and 
redshifts, infrared: galaxies 
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Introduction 



Less than a dec ade ago, observation s with the Submillimeter Common User Bolome- 
ter Array (SCUBA; iHolland et aD Il999h on the James Clerk Maxwell Telesco pe (JCMT) 



ident i fied a new and unexpected population of submillimeter galaxies (SMGs: iSmail et al. 



19971 ; iBarger et al.l Il998l ; iHughes et al.lll998l ). Subsequent surve ys with SCUBA and the 
Max-Planck Millimeter Bolometer (MAMBO; iKrevsa et aD Il998h array at the IRAM 30m 
telescope resolv ed a significant fra ction of the cosmic submillimeter background into individ- 



ual sources (see lBlain et al.ll2002l . and references therein). 



A detailed understanding of this population has emerged only slowly, due to their faint- 
ness at all short wavelengths and the difficult y of counterpart ide ntification. Photometric 
estimates of median redshifts around 2.5 — 3 (jCarilli fc Yunll2000l ) are consistent with the 
recent determination of a median redshift of z ~ 2 .2 for the < 50% of the population accessi- 



ble to optical spectroscopy (jChapman et al.ll2005l ). The optical redshifts of the radio /optical 
brigh t sub- c lass have been confirmed in ~ 15 ca s es through CO line d etections (jFrayer et al. 



19981 . Il99ft iNeri et al.l 120031 : iGreve et all l2005t iTacconi et all 120061 ). Despite all these ef- 



forts, the characterization of the redshift distribution remains incomplete, because of the 
large positional uncertainties for SMGs without interferometric counterparts and the large 
uncertainties in purely photometric redshift estimates. 

The tendency of SMGs to be faint in X-rays suggests that their larg e luminosities result 
primarily from high star formation rates ([Alexander et al.ll2003l . l2005bl ). Moreover, dynam- 
ical and gas pha se metallicity stud i es indicate that they have high (approaching 10 n M^) 



baryonic masses ( Genzel et al. 



Greve et al. 2005 



2003 



Tacconi et al 



Neri et al. 2003: Tecza et al. 2004: Swinbank et al. 2004 



20061 ). These high star formation rates and baryonic masses 
place SMGs at the assembly phase of massive galaxies at z = 2 — 3, in rapid and efficient 
starbursts likely fed by mergers. 

The contribution of AGN to the energy output of SMGs has strong implications for 
the origin of the cosmic submm background and the origin of th e correlation between black 



hole mass and spheroi d mass/velocity dispersion in local galaxies (IFerrarese &: Merritt 112000 



Gebhardt et al.l 120001 ) . Studies of the relationship between stellar and black hole mass in 
submillimeter galaxies, under the assumption of accretion at the Eddington rate, are con- 
sistent with a model where the supe rmassive black ho les in SMGs undergo rapid growth to 
reach the local M* — Mbh relation (IBorys et al.ll2005l ). 



SMGs thus mark a pivotal point in the evolution of massive galaxies and their central 
black holes. Observational characterisation of their redshift distribution, space densities, 
masses, metallicities, AGN content, and structure are all needed to understand their posi- 
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tion in the hierarchical growth of structure and the growth of massive galaxies by accre- 
tion of gas and merging of smaller galaxie s. Models of hierarchical galaxy formation (e.g. 
Kauffmann et al.l Il999l ; iBaugh et al.l 120031 ) are in the process of adapti ng to the new ob 



Baugh et all 12005 



servational constra ints including properties of the SMG population (e.g. 
Bower et al.ll2006l ). Extreme objects, such as the SMGs, can trace the formation of th e 
10 11 Mq galaxies already fully assembled at redshifts of z = 1.6 - 1.9 Jcimatti et all 120041 ). 
SMGs may also meet the constraint on rapid formation of low re dshift massive ellipti cals 
inferred from measurements of the a/Fe element abundance ratios ( Thomas et al.ll2005l ). 



Even at low resolution, rest-frame mi d-infrared spectra of gala xies can discriminate be- 



tween star formation and accretion (e.g. iGenzel fc Cesarskyl |2000| . and references therein) 



on the b asis of four distinc t spectral components observed in d usty galaxies in the local 

strong emission from 



Genzel et al. 


1998; 


Laurent et al. 


2000; 


Tran et al. 


2001) 



the polycyclic aromatic hydrocarbon (PAH) features found over a very wide range of star 
forming galaxies; (ii) a variable but usually small contribution of an HII region continuum 
steeply rising in the 6 — 15 //m rest wavelength range; (iii) a flatter PAH-free AGN con- 
tinuum, sometimes accompanied by addit ional A > 10 fim AGN related silicate emission 
(ISiebenmorgen et al.l 120051 ; lHao et al.l 120051 ) ; (iv) absorption featur es in the 6 — 8 /im range 
as well as the 9.7/im silicate absorption feature (jSpoon et al.ll2004l ). 



The anticorrelation between PAH feature strength relative to the AG N continuum and 
the ionization state of the ionized gas (IGenzel et al.lll998l ; iDale et al.ll2006l ) strongly supports 
these low resolution diagnostics. Using the Spitzer Space Telescope, the diagnostics based on 
the mid-infrared spectral components of dusty luminous and ultraluminous galaxies in the 
local universe can now be applied to dusty SMGs at high redshift. Our understanding of 
SMGs can be improved in three ways: 

Verification or determination of redshift. Mid-infrared spectral features, in par- 
ticular the narrow aromatic PAH features, if present, allow reasonably accurate redshift 
measurements (Az < 0.1) even for targets that are very faint at optical wavelengths. They 
can verify optical redshifts in cases where these are uncertain due to the faintness of the 
optical counterpart, due to th e presence of multiple candidate counterparts (e.g. 8 of 73 
SMGs in the redshift study of IChapman et al.l (120051 ) have multiple radio/optical counter- 
parts), or due to uncertain optical line identifications. Mid-infrared emission shares with 
CO line emission a reduced risk of an erroneous redshift assignment, because these tracers 
are more closely linked than the rest frame UV to the rest frame submm/far-infrared emis- 
sion that dominates the SMG's bolometric luminosity. They are thus less likely to measure 
the redshift of a misidentified neighbouring source. The difficulties of identifi cation with an 



optical/near-infrared source can be seen for example in the case of HDF850.1 (IDownes et al. 
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19991 ; iDunlop et al.ll2004l ). For sources without known optical redshifts, new mid- infrared 
spectroscopic redshifts strongly reduce the template ambiguities that are always possible 
for model-dependent photometric redshifts. The advantage of the mid-infrared spectra over 
CO is that, with existing instrumentation, there is a larger fractional bandwidth coverage 
than for mm spectroscopy. The disadvantage is that for a possible z > 4 tail of the SMG 
redshift distribution, the main PAH features leave the wavelength range of sensitive Spitzer 
spectroscopy. 

Evaluation of the relative importance of AGN and star formation using the rel- 
ative strength of PAH and continuum. Submillimeter sources are likely starburst-dominated 
but more l uminous than the threshold above which most local infrared galaxies are AGN 
dominated (jRigopoulou et al.lll999t IVeilleux et al.lll999l ; iTran et al.ll200ll ). Furthermore, al- 
though massive ellipticals in for mation should also be forming massive black holes in order 
to produce the Mrh — a relation (ITremaine et al.ll2002l ). evidence for energetically dominant 
AGN is scarce ( Alexander et al.l 120031 ) . Evaluating the ratio of starburst and AGN contri- 
butions will also allow us to investigate trends with other quantities. We adopt a ratio of 
7.7 fim PAH feature to local continuum of 1 as the border be tween starburst and AGN domi- 
nance in the bolometric luminosity, following the approach of lGenzel et al.l (119981 ) that is well 
matched to the rest wavelength coverage and S/N of our data. The rest-frame mid-infrared 
continuum, if isolated from the non-AGN components of the mid-infrared emis sion, can pro- 



vide constraints on the role of AGN, including highly X-ray absorbed sources (IKrabbe et al. 



200ll ; iLutz et al.l 120041 ) . Mid-IR spectra tracing the re-radiation of absorbed AGN emission 
thus also constrain the presence of Compton-thick AGN that are hard to separate from in- 
active objects in current high redshift X-ray data, because their X-ray photon statistics can 
be insufficient for clear identification of reflected AGN emission, or because the AGN may 
be fully covered. 

Constraining physical conditions in the starburst. By analogy with local galax- 
ies, dilute/cool regions are expected to show a larger PAH contribution to th e total infrared 



emiss ion and a cooler rest frame far-infrared peak than denser starbursts (e.g. iDale &: Helou 
20021 ). Very compact regions can show absorption features like those of water ice and hy- 
drocarbons. Some of the most luminous local ULIRGs exhibit strong con tinua with super- 



posed emission fe atures deviating from a canonical PAH feature shape (ITran et al.l 12001 
Spoon et al.l 120041 ) . Similar features may be found in SMGs. 



Spitzer spectra are also important to test the popular assumption (e.g. iHughes et al. 



19981 ) that the overall spectral energy distributions (SEDs) of SMGs are similar to those of 
local star forming Ultraluminous Infrared Galaxies (ULIRGs). 

In §2 we discuss the properties of our sample and in §3 observations and data analysis. 
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§4 presents the SMG spectra and discusses the results of their spectral decomposition. In §5 
we discuss the implications for the three main questions raised above. 



2. Sample Selection 



Among the several hundred SMGs detected up to now, our sample comprises exclusively 
bright (5*850 ^ m > 4.5 mJy) objects from blank field and cluster lens assisted surveys that 
have accurate positions from follow-up radio or mm interferometry. Objects from radio pre- 
selected surveys are not included because they may be biased in redshift and AGN content. 
No constraint on mid-infrared photometric flux, which could potentially introduce biases on 
SED and energy source, has been applied to this Cycle 1 sample. 

The 13 targets in our sample (Tab. [I]) cover the full range of properties of currently 
known bright SMGs. The ratio of optical/near- infrared and submm fluxes, for example, 



varies by almost two order s of magnitude (Tab. [3]). In our sample there are 6 ob 



published optical redshifts (IBarger et al.lll999t iLedlow et al.ll2002l ; I Chapman et al 



particular 3 cases where the redshift is confirmed by CO interferometry (jNeri et al 



ects with 



20051), in 



2003 



Greve et al.ll2005l ). and 7 sources for which no optical redshift is available due to their faint- 
ness in the optical/NIR. Including these objects in our study reduces biases that may affect 
a sample with measured optical redshifts only. The redshift distribution of a sample with 
optical redshifts will tend to avoid the optical "spectroscopic desert" (1.2 < z < 1.8). It 
may also favour sources with unobscured or mildly obscured AGNs because of their strong 
rest frame UV line emission. We include sources selected at 850 ^m (SCUBA) as well as 
sources selected at 1.2mm (MAMBO), to minimize selection effects due to redshift or dust 
temperature. The sample encomp asses the fi v e brig htest MAMBO sources with interfero- 
metric positions known before the IVoss et al.l (120061 ) observations, and a number of bright 
SCUBA sources with and without redshifts. 



3. Observations and Data Analysis 

We obtained low resolution (A/AA ~ 6O-12C0) long slit spectra using Spitzer IRS 
( IHouck et al.l 120041 ) in the staring mode. The detector is a 128 x 128 Si:Sb (LL1 and LL2 
modules) or Si:As (SL1 module) Blocked Impurity Band (BIB) array. 

The rest wavelength range required to detect mid-infrared PAH features and/or silicate 



AA is approximately constant as a function of A 
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absorption, encompasses at a minimum wavelengths from 5 to 10 /im. Observations are 
summarized in Tabled! Most of the sources were observed in the LL1 19.5 — 38.0 /mi module 
for 30 cycles of 120 s ramp duration. The total on-source integration time was 2 hours from 
addition of the independent spectra created by the telescope nods along the slit. For objects 
with z < 2.7 and most objects with unknown redshifts we added the LL2 14.0 — 21.3 /im 
range with 15 cycles of 120 s to ensure the required rest wavelength coverage for all plausible 
redshifts. In these cases, the slit on-source integration time for LL2 was 1 hour. Because 
of its low redshift, SMMJ16369+4057 was observed in the LL2 range for 10 cycles of 120 s, 
as well as in the SL1 7.4 — 14.5 /xm module for 5 cycles of 240 s (2400 s total on-source 
integration time). With this observation plan, we could in principle detect both the 6.2 /im 
and 7.7/tm PAH features of a starburst-like spectrum for a redshift range 1.4 < z < 3.8. For 
lower redshifts, we would observe only the 7.7 jum maximum, which fully exits the observed 
window at z ~ 1, but the longer wavelength 11.3 and 12.7 /im PAH features might still 
provide useful information in such a case. 

We reduced the data as follows. We subtracted, for each cycle, the two nod positions 
of the pipeline 14.0.0 basic calibrated data frames. In the difference just calculated, we 
replaced deviant pixels by values representative of their spectral neighborhoods. We sub- 
tracted residual wavelength dependent background, measured in source-free regions of the 
two dimensional difference spectra. In averaging all the cycles of the 2-dimensional sub- 
tracted frames, we excluded values more than three times the local noise away from the 
mean. The calibrated 1-dimensional spectra for the positive and the negative beams were 
extracted using the optimum extraction mode of the SPICE analysis package (version 1.4.1), 
and the two 1-dimensional spectra averaged in order to obtain the final spectrum. 

Since our measurements are close to the sensitivity limit of the instrument, for all data 
analysis we cut away the long and short wavelength ends of each module, where the noise is 
much higher and no signals are detected. 

In order to do a formal x 2 template fit for our sources, it is essential to determine the 
uncertainties on each point of the spectra. We used a processed and deglitched image of one 
of the undetected sources (the result was comparable for all of them). The uncertainty for 
each spectral point was calculated as the noise for each pixel times the square root of the 
number of pixels in the resolution element, which increases with wavelength. The noise for 
each pixel has been assumed to be the clipped standard deviation on a row of 20 pixels, all at 
the same wavelength. The resolution element is defined by the IRS pipeline, subdividing the 
region of the array covered by the spectrum into a sequence of trapezoid-shaped elements. 
This estimate may not fully reproduce the absolute noise level of a given observation, due to 
variations in strength of zodiacal light with position and epoch, but is a good approximation 
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of the change of noise as a function of wavelength. 



Results 



The spectra of the nine detected sources and the four non-detections are shown in Fig. [TJ 
As quantified by spectral fits below, the detections can be well reproduced by combinations 
of PAH features and continua. In one or two cases a heavily absorbed continuum is a possible 
alternative. We do not see evidence for strong continua like those seen i n local QSOs (which 



also s how silicate em i ssion outside our rest wavelength coverage, see ISiebenmorgen et al. 



2005; 


Hao et al. 


2005; 


Sturm et al. 


2005) 



ported by the PAH-dominated average of the SMG spectra (Fig. EI), which shows much larger 



PAH equivalent widths than those seen in QSOs (jSchweitzer et al.l 120061 ). 



In spectra of faint sources with poor S/N, only spectral features with sufficiently large 
equivalent widths can be used for redshift determination. At one extreme, these features 
are the strong PAH emission features of starburst galaxies. These features a ppear wheneve r 



the interstellar medium is exposed to moderately intense UV radiation (e.g. iDraine 1 120031 ) . 



The other extreme shows strong absorption features, the strongest being silicat e absorption. 



Exam p les of PAH featur e s in high- redshift ga laxi es observed with th e IRS are in lTeplitz et al. 



fl2005l ). lYan et al J (12005( 1. iLutz et al.l fl2005al) and 
local source is IRAS F00183-7111 JTran et al. 



Desai et al.l (120061 ) . An extremely absorbed 



20011 ). whose IRS spectrum is shown in 



Spoon et al.l (l2004f). Similar spec t ra have been observed in optically obscured high redshift 
24 jum sources ( jHouck et al.ll2005l ; IWeedman et al.l 120061 ). 



Redshifts can be determined by using either the set of strong PAH emission features 
or, for absorbed spectra, the 8 /im maximum and nearby silicate absorption. The strongest 
PAH feature is at 7.7/im (rest frame), so a similar redshift (Az ~ ±0.1) would be derived 
even if it is ambiguous whether the strongest feature is the 8 /im maximum or true PAH 
emission. The physical interpretation of the source, however, would be very different for the 
two alternatives. In order to correctly identify a feature as a PAH, we require an indication 
that the 6.2 /im PAH feature is present with the correct shape and (relative to the 7.7/im 
feature) flux. In our sample there is no case of ambiguity in identifying a feature with the 6.2 
or 7.7 /im PAH feature, since either both are detected or the large observed width indicates 
that it cannot be the narrow 6.2 /im PAH feature. For one of our sources, template fits 
identify such a single broad maximum with the 7.7 /tm PAH feature but, depending on rest 
wavelength range and S/N, it is also conceivable to identify such a single broad peak with 
the 8 /tm maximum of absorbed spectra. 
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We estimate redshifts by \ 2 fitting a selection of templates to the full spectra of the 
detected sources. Two of the t emplates are absorbed, likely AGN dominated, sources: NGC 



4418 and IRAS F00183 -7111 flSpoon et a 



200l|, |2004j) . As a third template we use the star- 



burst spectrum of M82 (jSturm et al.ll2000l ). both as observed (with the spectrum dominated 
by pure PAH emission) and with superposition of an additional continuum that is assumed to 
vary linearly over the short wavelength range covered. This last spectrum is similar to those 
of local universe infrared luminous galaxies having significant contributions to their bolo- 
metric luminositie s from both star formation and powerful AGN (e.g. NGC7469, Mrk273, 
Genzel et al.lll998l ). More complex fitting schemes have been used to interpret low resolution 
mid-infrared spectra, for e xample by allowin g additional obscuration to the PAH dominated 



starburst component (e.g. iTran et al.l l200ll ). or by decomposing the PAH component into 



individual feature s that can be well approximated by Lorentzians (e.g. lBoulanger et al.lll998 
Smith et al.ll2006l ). We do not adopt such schemes here because of the S/N of our data and 
their rest wavelength coverage, which is usually limited to shortward of the 9.6 /im silicate 
absorption by the IRS wavelength range. 

Results of the fits with all the templates are listed in Table El The best fit for each object 
is specified in Tab. Eland also shown in Fig. [H The redshifts derived from the best fit for each 
object are also listed in Tableland compared with known redshifts where available. IRS and 
previous redshifts are consistent for all the sources w ith previous redshifts from optical spec- 
troscopy (SMMJ09431+4700 bv lLedlow etldlhooi SMM.I 10519+5723. SMMJ10521+5719, 



SMMJ16369+4057, SMMJ16371+4053 bv IChapman et al. 



tional CO interferometry (SMMJ02 399-0136 by IFraver et al 



2005 ) an d sometimes from addi- 



19981. SMMJ 09431+4700 by 



Neri et al.ll2003l : iTacconi et al.ll2006l . SMMJ16371+4053 bv lGreve et al.ll2005[ ). Three sources 
had no accurately determined redshifts. For thes e, our new redsh i fts ar e inside the 68% con- 
fidence intervals of the photometric estimates of lAretxaga et al.l (120031 ) for at least three of 
the six different evolutionary models used. The formal fit uncertainties for the best fit are 
up to about Az = 0.02 but, because of the differences in derived redshifts depending on the 
template used (Tab. [2J), we assume an uncertainty Az = 0.06 for all the sources. For the 
six sources with both IRS and optical/CO redshifts measured we derive a reassuringly small 
standard deviation of only 0.014 for the difference between the best fitting IRS redshift and 
independent optical/CO redshift. For the uncertainty of a new IRS redshift from a spectrum 
similar in quality to our spectra, we nevertheless prefer the more conservative Az = 0.06 
which includes the possibility of template mismatch, since the smaller standard deviation for 
the six sources includes only spectra reasonably well fit by the M82 (plus continuum) tem- 
plates. Az = 0.06 is an overestimate if the identification of PAHs in the spectrum is beyond 
any doubt. Table [3] also classifies sources according to how their spectra are characterized: 
PAH emission (M82 best fit), silicate absorption or 8 fim maximum (IRAS F00183— 7111 or 



- 9- 



NGC4418 best fit) or a superposition of PAH emission and a linearly rising continuum. 

In the remainder of this section, we briefly describe the properties of the nine SMGs with 
new mid-infrared rest-frame spectra. We also discuss the four non-detections and constraints 
on their possible redshifts and/or SED properties. 

SMMJ00266+1708 The best fit for this source is the combination of M82 PAH tem- 
plate and weak AGN continuum (see Fig. [T]), but the \ 2 test gives similar results also for 
the absorbed continuum NGC 4418 template (see Tab. [2]). The spectrum lacks an evident 
6.2 /im PAH feature, the wavelength of which still falls in the observed band. We proceed 
with the best fitting PAH plus continuum interpretation but note that a classification of this 
source as an absorbed A GN is clearly no t exclu ded. The previous photometric redshift of 
z = 2.7±oi estim ated by lAretxaga et al.l (120031 ) is consistent with our spectroscopic value 



2.73 ± 0.06. iFraver et al 



(120001 ) identify a faint ERO counterpart for this object and 



estimate a lensing magnification 2.4 ± 0.5, which we adopt for our analysis. 

SMMJ02399— 0136 The mid- infrared spectrum of this source was already presented 
by lLutz et al.l (j2005al ). It contains well detected 6.2 /im and 7.7/xm PAH features superposed 
on a strong continuum (see Fig. [1]). This source is at the transition between predominantly 
starbu rst powered a n d pre domi nantly AGN powered, according to the mid-infrared diagnos- 



tics of Genzel et a 



optical spectrum (llyison et al. 



X-ray observations (IBautz et al 



(1998 1 and lLaurent et al.l (120001 ). AGN signatures are also seen in the 



998). Combining mid-infrared spectroscopy and Chandra 



2000l ) , we conclude that this source is powered by rou ghly 



equal contributions of star formation and a Compton-thick AGN (ILutz et al.l l2005aj ) . in 
agreement with constraints on the importance of star formation in this object from molec 



ular gas mass and its position relative to the radio/far-infrared correlation (IFraver et al. 



19981 ). Its lu minosity (see Tab. EH) is calculated taking into account a lensing magnification 
factor of 2.5 Jlvison et al.lll998h . 



SMMJ09429+4659 This source is well fitted with the starburst-like M82 spectrum. 
No accurate redshift was pr eviously known for th is object. Using the radio-sub mm spectral 
index and the models from ICarilli fc Yunl (120001 ). the estimated redshift is z = 0.4 ± 0.3. 
The intense radio emission suggests, though, that there is a significant radio contribution 
from an AGN, even if we do not see AGN emission in the mid-infrared spectrum, making 
the value estimated from the radio-submm relation only a l ower limit. Assum ing this galaxy 
follows the K — z relationship fo r powerful radio galaxies (jJarvis et al.l 120011 ). it most likely 
lies at z > 2 (ILedlow et al.l 120021 ) . Our spectroscopic value of z = 2.38 ± 0.06 confirms the 
high redshift of this object. The luminosity of this source (see Tab. E2) is calculated by taking 
into account a lensing magnification factor of 1.3 ( Cowie et al. 20021 ). 
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SMMJ09431+470 We pointed I RS a t the mm position flNeri et all 120031 ) of compo- 
nent H7 in the notation of lLedlow et all (120021 ). but note that component H6 is included in the 
observing aperture as well. This source has the highest redshift of the sample (z = 3.36), not 
accounting for the non-detected objects (see discussion below). It is well fitted with a PAH 
spectrum plus a strong linearly rising continuum, so we infer that it is powered by both star- 
burst activity and a power ful AGN (see Fig. [D) . The optical spectrum of component H6 shows 
features of a weak AGN (ILedlow et al.ll2002l ). The spectral properties, line widths and line 
ratios of this galaxy a re very similar to those seen for narrow- line Seyfert 1 galaxies (NLSyls; 
Crenshaw et aDll99lh . XMM-Newton observations (ILedlow et al.ll2002l ) suggest that the in- 



trinsic X-ray luminosity of the AGN is modest (L 2 _ 



10kcV 



< 10 44 erg s 1 cm 2 ) 



unless there is 



heavy obscuration. Comparison to the strong rest frame mid-IR continuum seen in the IRS 
spectrum suggests the latter is indeed the case: at L 2 -iokcv/ zy L v (6 fim) < 0.015, this source 
falls more than an order of magnitude below t he relation betwe en unobscured 2 — lOkeV 
luminosity and 6 fjm continuum for local AGN jLutz et al.l l2004l ) . SMMJ09431+4700 hosts 
a heavily obscured or Compton-thick AGN, the location of which we cannot firmly ascribe 
to component H6 or H7. The luminosity of SMMJ094 31+4700 (see Tab . [3]) is calculated 
taking into account a lensing amplification factor of 1.2 (jCowie et al.ll2002l ). 



SMMJ10519+5723 The lowest \ 2 fit for this object is the one with the M82 starburst 
template, which is supported by the tentative detection of a 6.2 /im PAH feature. Still, the 
X 2 is not much higher for the absorbed continuum template NGC 4418. We interpret this 



sourc e as powered by starburst activity, consistent with optical spectroscopy (jChapman et al. 



20051 ). but note that an absorbed continuum interp retation cannot b e firml y excluded. The 
spectr oscopic redshi f t sugg ested for this source bylChapman et al.l (j2003al ) was z = 3.699. 
Later, Egami et al. ( 2004 ) and Chapman et al. ( 20051 ) indicated values of z = 2.69 and 
z = 2.686 respectively. Our fits give z = 2.67 ± 0.06 using the starburst template and 
z = 2.64 ± 0.06 using the obscured one. Both redshifts are in agreement with the latest 
optical results. 

SMMJ10521+5719 The spectrum of this source shows PAH features plus a weak con 



tinuu m. The source does not contain AGN features in the optical spectrum (jChapman et al. 



20051 ). The mid- infrared AGN continuum is detected but the feature-to-continuum ratio is 
much higher than in the cases of SMMJ02399-0136 and SMMJ09431+4700, thus suggesting 
a smaller AGN contribution. The presence of an AGN is confirmed as well from radio a nd 
X-ray emission that show the possible presence of a radio loud quasar (jlvison et al.ll2002l ). 



MMJ154127+6 616 The mid-infrared spectrum of this source was already published 
by iLutz et al.l (j2005al ). We reprocessed the data using a later version of the IRS pipeline 
(14.0.0). The new reduction further increases the similarity to a starburst spectrum with 
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well defined 6.2 /mi and 7.7/mi PAH features and is well fitted by a M82 spectrum plus 
a very weak flat continuum. The y 2 is not much worse adopting the NGC4418 template 
which better matches the emission at long wavelengths, but the clear presence of both PAH 
features strongly increases confidence that this source is powered by star formation. 

SMMJ16369+4057 This source has the lowest redshift of the sample (z = 1.21) and is 
w ell fitted with the M82 spectrum plus a very weak flat continuum. The optical spectroscopy 
of IChapman et al.l (120051 ) detects typical starburst lines, consistent with our result. 



SMMJ16371+4053 The mid-infrared spectrum of this source shows clear 6.2 /mi and 
7. 7 /mi PAH features. A flat continuum is detected. The mid-infrared spectrum shows that 
it is powered mainly by starburst activity, even though an AGN is probably a lso present. 
AGN lines have been detected in the optical spectrum of I Chapman et al.l (120051 ). 



IRS non-detections The sources for which we have neither detected features nor con- 
tinua in the IRS spectra despite accurate interferometric positions are MMJ120517— 0743.1, 
MMJ120539-0745.4, MMJ120546-0741.5 and MMJ154127+6615. Multiple arguments sup- 



port the reality of these mm sources: they are wel 



(Ber toldi et al.ll2000l ; iDannerbauer et al 



2002 



20031 ) . have weak VLA counterparts (jBertoldi et al.1 



detected in the original M AMBO data 



2004 ), are confirmed bv SCUBA (|Eales et al. 



2000; IDannerbauer et al.ll2004) and have 



dust continuum emission direct ly confirmed and located by mm interferometry (IDannerbauer et al 
20021 . |2004| ; lDannerbauerll2004l ) . The faintness or non-detection of their optical/near infrared 
counterparts (see also Table 3) lead these authors to the conclusion that these sources must 
be at very high redshift and/or highly obscured. 

We now consider each of these scenarios in turn. If the IRS non-detections are intrinsi- 
cally similar to the detected sources but simply more distant, we can place one lower limit 
on their redshifts by requiring that the 7.7 — 8.6 /im PAH complex or the 8 /mi maxima in 
any absorbed spectra have remained undetected because they were shifted to wavelengths 
longer than ~ 35 /mi, where the detector sensitivity falls rapidly and the noise increases. 
This argument implies redshifts greater than 3.6. Similarly, since 6.2 /mi features in PAH 
emission spectra would already be lost in the noise at wavelengths longer than ~ 30 /mi (be- 
ing somewhat weaker), such sources could only be IRS non-detections for z > 3.8. Hence, 
we would conclude for all undetected sources redshifts > 3.6. 

If the IRS non-detections lie at redshifts similar to those of the IRS detections, their 
SEDs must be significantly different. We infer in Section 5.2 below that the detected sources 
already have a ratio of mid-infrared PAH to far-infrared emission similar to that of local 
ULIRGs. This implies that to escape IRS detection at the same redshifts, the non-detections 
would need to be even more heavily obscured in the mid-infrared than local ULIRGs. By way 
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of example, c onsider Arp 220, wel l known for its ex tremely low ratio of mid- to far- infrared 



Sanders et"aDll988l : Irlaas et al.lboOll ). At z ~ 2.5, an extreme Arp 220-like 



emission (e.g. 

SED scaled to the millimeter fluxes of the IRS non-detections would still manifest a broad 
PAH feature peaking at the 0.15 — 0.2 mJy level. No such feature is indicated in the spectra 
of the four non-detections (Fig. [p. 

Previous assessments of the "higher redshift" and "higher obscuration" scenarios for 
these sources have been based on faint or undetected near-IR counterparts. For MMJ120517- 
0743.1, MMJ120539-0745.4, and MMJ120546-0741.5, faint i^-band magnitudes and the 
assumption of SEDs sim ilar to those of local ULIRGs imply very high redshifts z > 4 
( IDannerbauer et al.ll2002l ). Similarly, the K > 21.2 counterpart of MMJ154127+6615 implies 
a redshift z > 3. This particular argument is weakened by evidence that SMGs can have rest- 
frame UV/optical obscurations greater than those of local ULIRGs; SMMJ00266+1708, for 
example, has a very faint near-IR counterpart but (based on our data) z = 2.73. However, if 
we generalize the argument to the rest-frame mid-infrared, which is more difficult to obscure 
than the UV/optical, we are on stronger ground. Tellingly, our IRS non-detections are at 
least ~ 3 times fainter in the rest-frame mid-inf rared than even SMMJ0026 6+1708. This 
result is confirmed by the mid-infrared imaging of ICharmandaris et al.l (120041 ). who place all 
four sources (three not detected and one tentatively detected) at the low end of the mid-IR 
to submm flux ratio distribution. 

A further reason to prefer the "higher redshift" scenario for the IRS no n-detections is 



that t he radio counterparts to a ll four are faint, despite their bright mm fluxes (IBertoldi et al. 



2000l ; IDannerbauer et al.ll2004j ). The me an of their ratio o f 850 ii m and 1.4 GHz flux densities 
(212 ± 77) is about twice that of the I Chapman et al.l (120051 ) sample (95 ± 9), and still 
higher than both the value of 108 ± 9, obtained for the same sample after exclusion of 
potentially radio-loud (Si.4GHz > 200 /iJy) SMGs, and the value of 150 ± 20, derived from 
a 15 sources subsample matching the mean SCUBA flux o f our four non-detections. This 
trend is consistent with higher redshifts (jCarilli fc Yunll2000l ). although the scatter about the 
radio/submm vs. z relation is large. Taken together, the various lines of evidence suggest 
that the four IRS non-detections do lie at high redshift, although with all arguments still 
based on SED assumptions. Direct spectroscopic redshifts will be needed for a definitive 
conclusion. 
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Discussion 



5.1. The redshift distribution has a median of z ~ 2.8 



Figure [2] shows the redshift distribution for our sample. Taking into account the lower 
limits adopted for the undetected sources, we derive a median redshift of z — 2.79 for the 
full sample of 13 SMGs. The median redshift of the 9 detected sources is z = 2.69. Th ese 
values are noticeably higher than the median z = 2.2 measured by I Chapman et al.l (120051 ) for 
their sample of 73 submillimeter galaxies with optical redshifts, and also higher than their 
estimate of z = 2.3 for the extrapolation to the full SMG population. None of our new IRS 
redshifts is in the optical "redshift gap " (1.2 < z 5, 1-8) of th e optical spectroscopic census of 



the submillimeter galaxy population (j Chapman et al. 



20051 ). where some new IRS redshifts 



might have been expected. It is natural to assume that the optical reds hifts are still biased 



towa rds the optically bright and low redshift part of the population (Da nnerbauer 



2004 ) . The recent discovery of a submillimeter galaxy at redshift z ~ 4 fcnudsen et alJbood ) 



et al. 



si milarly indicates an exte nsion of the S MG redshift d istrib ution beyond the one established 
by IChapman et al.l (120051 ). In contrast, iPope et al.l (120061 ) suggest from a combination of 
spectroscopic and photometric redshifts a median redshift of 2.0 for those SMGs in the HDF- 
North region that they consider securely identified through radio or Spitzer counterparts. 



Our results suggest a modest extension of the IChapman et al.l (120051 ) redshift distribu- 
tion towards a larger high redshift tail, but the number statistics of our and other current 
SMG samples with redshifts is small. To investigate the significance of this difference, we 
have run simple Monte-Carlo simulations to estimate the probability of obt aining by chance 



Chapman et al. 



Chapman et al. 



a med ian redshift > 2.79 for a 13 objects sample drawn randomly from a 
(2005) redshift distribution. Using the distributions plotted in Fig. 4 of 
(120051 ). this probability is a low 0.3% drawing 13-objects samples from their optical SMG 
redshift distribution with overall median 2.2, and a still interesting 8% for drawing from 
their suggested extrapolation to the overall SMG redshift distribution with median 2.3. The 
latter comparison is conservative in ignoring the fact that our sample with mostly radio- 
interferometric identifications is still biased against the very highest redshift objects. 

In addition to these simple statistical comparisons, and perhaps more important, there 
can be effects of field-to-field varia tions including sp ikes in the redshift distribution for the 
current small area SMG surveys (jBlain et aL 20041). A potent i al z > 4 redshift spike in 



the NDF region observed with MAMBO by iDannerbauer et al.l (120041 ) from which three of 
our IRS non-detections are drawn, for ex ample , could emph asize the high z component in 
our small sample. Likewise, the results of IPope et al.l (120061 ) could have been driven in the 
opposite direction by being drawn from a single field. These constraints clearly call for 
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further analyses using larger samples from large and widely separated fields. 



5.2. SMGs have ULIRG-like SEDs and are largely starburst-powered 



The rest frame mid-infrared spectra of SMGs and their comparison to the far-infrared 
part of the SEDs hold interesting clues about their physical properties, energy sources and 
radiation fields. The average spectrum of the nine detected SMGs, individually scaled to 
the same rest wavelength flux S , 222^m = 15mJy to give all sources equal weight (see Fig. [3]), 
clearly shows the PAH features at 6.2 and 7.7 /im but relatively weak continuum. The 
rest wavelength of 222 /im was chosen for this scaling because it is well constrained by 
observations, since it is the rest frame wavelength for the SCUBA 850 /iin flux at redshift 
z = 2.8, about the median redshift of the sample. Scaling by the rest frame far-infrared 
emission is close to scaling by bolometric flux, which is appropriate for interpreting the mid- 
infrared diagnostics in Fig. [3] in terms of energy sources of a typical SMG without biasing 
towards star formation (PAH) or AGN (continuum). For galaxies at different redshifts, the 
rest frame 222 /im continuum emission was extrapolated from the SCUBA flux. For this 
step as well as in quantifying the bulk properties of our targe ts, we assign luminosities and 
temperat ures following the approach of I Chapman et al.l (120051 ) who use the local FIR/radio 
relation ( Helou et al.lll985l ) to assign FIR SEDs to sources with known redshift, radio, and 
submm fluxes. The adopted far-infrared SED model is a grey body with S U) t oc k(v)B V)T 
with k(u) oc v 13 and (3 = 1.5. We calculated Td fo r our galaxies us i ng two photometric points 
(850 /im and 1.4 GHz) and the relation found by lChapman et al.l (120051 ) 

1 + z 



T d oc 



($850 Mm/ 5l.4GHz) 



0.26 



1 



with a proportionality constant of 6.29 (see Tab. [3]). Infrared luminosities are calculated 
as the integral between 8 and 1000 /mi of the SED, assuming a ACDM cosmology with 
Hq = 70 kms _1 Mpc _1 , fl,M = 0.3 and Q\ = 0.7 (see Tab. [3]). Like other SMGs, our sample 
objects are inferred to be very luminous objects (L IR ~ 10 13 Lq). 

The spectrum of Fig. [3] pr ovides a first and direct indicati on that our sources are, on 
average, starburst-like (see also iMenendez-Delmestre et al.ll2007l . for a similar conclusion for 
several lower redshift SMGs). It can be seen as a superposition of a M82-type PAH spectrum 
and a weak additional continuum. Comparison with the spectra of t he mostly starburst- 
dominated local ULIRG population (see e.g. Fig. 1 in lLutz et al.lll998l ) suggests that SMGs 
are scaled up versions of these objects. A further proof of this comes from a comparison of 
the ratio of PAH features and S222fim continua. Figure H] shows the ratio of peak flux density 
of the 7.7 /im PAH feature, after continuum subtraction, to the continuum flux density at 



- 15 - 



222 /im rest frame for all the detected sources of our sample and for 11 local ULIRGs with 
PAH emission a nd good FIR photometry . The ULIRG data have been taken from ISOPHOT- 
S observations (IRigopoulou et al.l I1999T ) and the continua from slight extrapolations of the 
far- infrared photometry of iKlaas et al.l (120011 ). which extends to an observed wavelength of 
200 /im. The PAH-to-far-infrared ratios of our sources are fully consistent with those of 
the local ULIRG population. These SED properties are in agreement with the conclusion, 
from spatially resolved mm interferometry, that SMGs are similar to local ULIRGs suitably 
scaled for the ir larger masses, lum inosities and star formation rates, as well as their greater 
gas fractions (jTacconi et al.ll2006l ). 



Given this similarity of SMGs and local ULIRGs in the comparison between PAH and 
far- infrared parts of the SED, it is instructiv e to also compar e the dust temperatures for 
z ~ 2.5 SMGs with those of local galaxies. I Chapman et al.l (120051 ) infer Td ~ 36 ± 7K 
for SMGs which can be compared to local galaxies in two ways: (1) SMGs have colder 
dust temperatures than local ULIRGs of t he same luminos it y (tha t means colder dust than 
so-called HyLIRGs). For local HyLIRGs, Ichapman et~aD J2003bh find T d ~ 42K at L ~ 
10 13 Lq. This could be due t o their difference in energy sources: local HyLIRGs seem 
predominantly AGN driven (e.g. lTran et al.ll200ll ). while many SMGs of the same luminosity 
are starbursts. (2) SMGs have similar dust temperatures to the bulk of t he local ULIRG 



popula tion which has luminosities just above 10 12 Lq. For local ULIRGs, iChapman et al. 



(j2003al ) infer Td ~ 36K. This comparison is in line with other indications that SMGs are 
"scaled up" versions of such ULIRGs rather than direct analogs of local HyLIRGs. 

From the rest frame mid-infrared spectra it is possible to evaluate the presence and the 
strength of a possible AGN contribution to the very large IR luminosity of the SMGs of our 



samp 



2003 



e. Ma ny SMGs show evidence of an AGN, fr om X-ray observations ([Alexander et al. 



2005a 



bf) or opti cal/near-infrared spectra (e.g. llvison et al.lll998l ; ISwinbank et al 



2004 



Chapman et al.l 120051 ) . However, the most important question is not whether there are 
detectable AGN signatures, but whether or not the AGN is a significant contributor to the 
luminosity of the galaxies. From our spectral decompositions, we can use the ratio of PAH 
7.7 /im peak and local continuum as an indicator of the AGN and star for mation contribution 



to the bolometric luminosity. Followin g studies of local infrared galaxies (IGenzel et al.lll998 



Laurent et al.l l2000l : iTran et al.l l200ll ). we adopt a feature to continuum ratio of 1 as the 
border between predominantly star formation and predominantly AGN powered^]. The fitted 



2 Our decomposition by template fit differs from the one used by Genzel et al. ( 19981) . which interpolates 
between two points in the observed spectrum to define the continuum, in approaching very high feature to 
continuum values in the limit of almost pure star formation. The two methods give very similar results in 
the presence of significant AGN continuum 
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fluxes for the 7.7 /zm PAH peak and local continuum are listed in Table [3) 



In our sample, there is no trace of an AGN continuum in the adopted fits to the rest frame 
mid-infrared spectrum for two sources (SMMJ09429+4659 and SMMJ10519+5723). Seven 
sources (SMMJ00266+1708, SMMJ02399-0136, SMMJ09431+4700, SMMJ10521+5719, 
MMJ154127+6616, SMMJ16369+4057, SMMJ16371+4053) are well fitted by a superposi- 
tion of an AGN continuum and PAH features. We note again that the alternative fit by an 
obscured (AGN?) continuum cannot be firmly excluded for one of these (SMMJ00266+1708). 
Assuming that a feature to continuum ratio of 1 means similar bolometric contribution from 
star formation and accretion, we conclude that, of these nine targets, two are pure starbursts 
with at best very weak AGN, five have AGN with modest contributions of the order 20%, and 
two (SMMJ02399-0136, SMMJ09431+4700) have strong AGN contributing slightly above 
half of the bolometric luminosity. We cannot constrain the starburst or AGN nature of 
the four undetected sources that are likely at high redshift, apart from stating that their 
non-detection means absence of a strong unobscured hot AGN continuum even for redshifts 
somewhat above 3.6. Our sample is small, so it seems reasonable to expect that the full sub- 
millimeter population will show an even wider range of AGN pro pertie s perhaps in c luding 
less obscured, energetically dominant AGNs. lEgami et al.l (120041 ) and llvison et al.l (120041 ) 
have used mid-IR photometry as a diagnostic tool to put limits on the AGN contribution 
to the IR luminosity and agree that fewer than 25 % of the SCUBA/MAM BO sources ob- 
served are AGN powered. From X-ray observations. Alexander et al.l (l2005bl ) found that, on 
average, the AGN contribution in a large sample of radio-detected SMGs was likely to be 
modest (~ 10%). This result assumes that SMGs do not have a substantially larger dust- 
covering factor than optically selected quasars, and that there is no significant number of 
fully Compton-thick AGNs that are hard to detect in X-rays. Our finding of starbursts be- 
ing prevalent in the mid-infrared spectra extends this result by showing that SMGs typically 
do not contain such dominant X-ray obscured AGNs with strong mid-infrared continuum 



re-emission. Obscurec 
infrared surveys ( e.g. 



AGNs with strong mid-infrar ed continuum are found in Spitzer mid- 



Martmez-Sansigre et al.ll2005l ). but show only little overlap with the 



SMG population (ILutz et al.ll2005bl ). In general it seems that the AGN contribution to the 
IR luminosity of most SMGs is small when compared to heating from star-formation activity. 

While the IRS spectra argue against dominant AGN being typical for SMGs, they 
can still help to constrain the obscuration of the lesser AGN found, independent of whethe r 
sufficient photons are available for a detailed X-ray spectral analysis ([Alexander et al.ll2005bl ) . 
From comparison to X-ray data, we have argued in the discussion of the individual sources 
for high X-ray obscuration of the two strongest AGN in our sample. For z = 2.8, a rest frame 
6 /im continuum of ~ 0.1 mJy for the weaker AGN, and the relation of rest frame 2 — 10 keV 
flux and 6 /im AGN continuum from lLutz et al.l (120041 ). a rest frame 2 — lOkeV emission of 
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4 x 10~ 15 erg s _1 cm -2 is expected for an unobscured AGN following this relation. This is in 
the sensitivity regime of current X-ray data in the corresponding observed frame band, and 
can constrain the obscuration of such "minor" AGN with the caveat that large samples are 
needed given the variations in intrinsic AGN SEDs and the corresponding significant scatter 
of the mid-IR/X-ray relation. We have used Chandra archival X-ray data to put limits on the 
observed frame 0.5 - 2keV emission of SMMJ00266+1708 (< 0.56 x 10" 15 erg s^cm" 2 ) and 
MMJ154127+6616 (< 0.16 x 10~ 15 ). From XMM-Newton observations of SMMJ10521+5719 
in the same soft band, a value of 0.16 x 10~ 15 erg s _1 cm -2 is derived (Brunner et al., in 
prep.). These limits and fluxes imply that the rest frame 2 — lOkeV emission is lower than 
the extrapolation from the mid-infrared continuum to unobscured X-rays. Noticeable X-ray 
obscuration may thus b e found in many of the minor AGN in SMGs, in agreement with 
Alexander etaD J2005bh . 



Metallicity and dust-to-gas ratios in many high redshift galaxies are expected to be lower 
than at low redshift. Low me tallicity systems show weaker mid-infrared PAH emission bands 
(e.g. lEngelbracht et al.ll2005l ). In addition to differences in radiation fields, this is probably 
due to the fact that these galaxies are young and thus may lack the carbon-rich AGB stars 
required to form the PAH molecules. However, the enrichment will proceed once intense 
star formation activity is underway for a sufficient time, or if there has been preceding star 
formation. This seems to be the case for SMGs. They show very clear PAH features as 
tracers of their intense starbursts. We can therefore assume that high metallicity is typical 
for the population of massive SM Gs, consistent with the supersolar metallicity derived from 
nebular emission fo r SMM14011 (ITecza et al.l 12004 ) and the roughly solar abundances in 
the SMG sample of lSwinbank et al.l tooi ). In fact, because of the intense star formation, 
the metallicity of these systems should r apidly approach that of their likely present-day 
descendants: luminous elliptical galaxies (ISwinbank et al.l 120041 ) . These observations and 
the scenario that SMGs evolve to ellipticals are in full agreement with the fossil record that 
the formation of the stars of local ellipticals must have happened rapidly and at high redshift 
jThomas et al.ll2005h . 



6. Conclusions 

We have presented Spitzer mid-infrared spectra of a sample of 13 submillimeter galaxies. 
For nine of them, we have unambiguous detections of PAH spectral features and/or mid- 
infrared continua that allow us to constrain energy sources in these objects and to determine, 
in three cases for the first time, their redshifts. 

The IRS detections alone have a median z ~ 2.7. If the four IRS non-detections lie 
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at similar redshifts, their rest frame mid-infrared obscurations would have to be even more 
extreme than those of local ULIRGs. More plausibly, the four IRS non-detections lie at 
higher redshifts (> 3.6), giving a median z ~ 2.8 for the full set of 13. Although our 
sample is small, this result may indicate an extension to higher redshift of the SMG redshift 
distribution relative to radio-preselected samples with optical redshifts. 

In the majority of cases, the detection of PAH emission and the weakness of AGN 
continua indicate that these galaxies are mainly starburst-powered. This result agrees with 
previous X-ray, optical and SED studies that indicate only a small AGN contribution to the 
IR luminosity compared to heating from star-formation activity. Our work extends these 
studies by also constraining the role of highly obscured AGN. 

The SED properties of our galaxies are in agreement with the SMGs being scaled-up 
versions of the compact star-forming regions in local ULIRGs. 

The existence of star formation dominated systems at infrared luminosities in excess of 
10 13 Lq is unique to the high redshift universe. The presence of high luminosity starbursts in 
SMGs may be related to their higher gas fractions (IGreve et al.ll2005l ; iTacconi et al.l 120061 ) . 



Mid-infrared spectroscopy with IRS, together with ancillary observations from the op- 
tical through radio wavelengths, can play a central role in understanding the nature of 
submillimeter galaxies and can be a powerful tool for probing the earliest and most dramatic 
stage of the evolution of galaxies. 
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following people for helpful discussions and support during the work: M.Righi, F.Braglia, 
L.Conversi. 



REFERENCES 

Alexander, D.M., et al. 2003, AJ, 125, 383 



Alexander, D.M., Smail, I., Bauer, F.E., Chapman, S.C., Blain, A.W., Brandt, W.N., Ivison, 
R.J. 2005, Nature, 434, 738 



-19- 

Alexander, D.M., Bauer, F.E., Chapman, S.C., Smail, I., Blain, A.W., Brandt, W.N., Ivison, 
R.J. 2005, ApJ, 632, 736 

Aretxaga, I., Hughes, D.H., Chapin, E.L., Gaztanaga, E., Dunlop, J.S., Ivison, R.J. 2003, 
MNRAS, 342, 759 

Barger, A.J., Cowie, L.L., Sanders, D.B., Fulton, E., Taniguchi, Y., Sato, Y., Kawara, K., 
Okuda, H. 1998, Nature, 394, 248 

Barger, A.J., Cowie, L.L., Smail, I.R., Ivison, R.J., Blain, A.W., Kneib, J.-P. 1999, AJ, 117, 
2656 

Baugh, CM., Benson, A. J., Cole, S., Frenk, C.S., Lacey, C.G. 2003, |astro-ph/0203051| 

Baugh, CM., Lacey, C.G, Frenk, C.S., Granato, G.L., Silva, L., Bressan, A., Benson, A. J., 
Cole, S. 2005, MNRAS, 356, 1191 

Bautz, M.W., Malm, M.R., BaganofT, F.K., Ricker, G.R., Canizares, C.R., Brandt, W.N., 
Hornschemeier, A.E., Garmire, CP. 2000, ApJ, 543, L119 

Bertoldi, F., et al. 2000, A&A, 360, 92 

Blain, A.W., Smail, I., Ivison, R.J., Kneib, J.-P, Frayer, D.T., 2002, PhR, 369, 111 

Blain, A.W., Chapman, S.C, Smail, I., Ivison, R. 2004, ApJ, 611, 725 

Borys, C, Smail, I., Chapman, S.C, Blain, A.W., Alexander, D.M., Ivison, R.J. 2005, ApJ, 
635, 853 

Bower, R.C, Benson, A.J., Malbon, R., Helly, J.C, Frenk, C.S., Baugh, CM., Cole, S., 
Lacey, C.G. 2006, MNRAS, 370, 645 

Boulanger, F., Boisssel, P., Cesarsky, D., Ryter, C. 1998, A&A, 339, 194 

Carilli, C.L., Yun, M.S. 2000, ApJ, 530, 618 

Chapman, S.C, Blain, A.W., Ivison, R.J., Smail, I.R. 2003, Nature, 422, 695 
Chapman, S.C, Helou, C, Lewis, G.F., Dale, D.A. 2003, ApJ, 588, 186 
Chapman, S.C, Blain, A.W., Smail, I., Ivison, R.J. 2005, ApJ, 622, 772 
Charmandaris V., et al. 2004, ApJS, 154, 142 
Cimatti A., et al. 2004, Nature, 430, 184 



-20- 

Colbert, J.W., et al. 1999, ApJ, 511, 721 

Cowie, L.L., Barger, A., Kneib, J.-P. 2002, AJ, 123, 2197 

Crenshaw D.M., Peterson, B.M., Korista, K.T., Wagner, R.M., Aufdenberg, J.P 1991, AJ, 
101, 1202 

Dale, D., Helou, G. 2002, ApJ, 576, 159 
Dale, D., et al. 2006, ApJ, 646, 161 

Dannerbauer, H., Lehnert, M.D., Lutz, D., Tacconi, L., Bertoldi, F., Carilli, C, Genzel, R., 
Menten, K. 2002, ApJ, 573, 473 

Dannerbauer, H., Lehnert, M.D., Lutz, D., Tacconi, L., Bertoldi, F., Carilli, C, Genzel, R., 
Menten, K. 2004, ApJ, 606, 664 

Dannerbauer, H. 2004, PhD Thesis, Ludwig-Maximilians-Universitat Miinchen 

Desai, V., et al. 2006, ApJ, 641, 133 

Downes, D., et al. 1999, A&A, 347, 809 

Drainc, B.T. 2003, ARA&A, 41, 241 

Dunlop, J.S., et al. 2004, MNRAS, 350, 769 

Eales, S., Bertoldi, F., Ivison, R., Carilli, C, Dunne, L., Owen, F. 2003, MNRAS, 344, 169 
Egami, E., et al. 2004, ApJS, 154, 130 

Engelbracht, C.W., Gordon, K.D., Rieke, G.H., Werner M.W.,Dale, D.A., Latter, W.B. 
2005, ApJ, 628, L29 

Farrah, D., Afonso, J., Efstathiou, A., Rowan-Robinson, M., Fox, M., Clements, D. 2003, 
MNRAS, 343, 585 

Ferrarese, L., Merritt, D. 2000, ApJ, 549, 215 

Forster Schreiber, N., Sauvage, M., Charmandaris, V., Laurent, O., Gallais, P., Mirabel, 
I.F., Vigroux, L. 2003, A&A, 399, 833 

Frayer, D.T., Ivison, R.J., Scoville, N.Z., Yun, M., Evans, A.S., Smail, I., Blain, A.W., 
Kneib, J.-P. 1998, ApJ, 506, L7 



Frayer, D.T., et al. 1999, ApJ, 514, L13 

Frayer, D.T., Smail, I., Ivison, R.J., Scoville, N.Z. 2000, AJ, 120, 1668 

Fukugita, M., Hogan, C.J., Peebles, P.J.E. 1998, ApJ, 503, 518 

Gebhardt, K., et al. 2000, ApJ, 539, L13 

Genzel, R., et al. 1998, ApJ, 498, 579 

Genzel, R., Cesarsky, C. 2000, ARA&A, 38, 761 

Genzel, R., Baker, A.J., Tacconi, L.J., Lutz, D., Cox, P., Guilloteau, S., Omont, A. 2003, 
ApJ, 584, 633 

Greve, T.R., et al. 2005, MNRAS, 359, 1165 

Haas, M., Klaas, U., Miiller, S.A.H., Chini, R., Coulson, I. 2001, A&A, 367, L9 
Hao, L., et al. 2005, ApJ, 625, L75 

Helou, G., Soifer, T., Rowan-Robinson, M., 1985, ApJ, 298, L7 PASP, 116, 975 
Holland, W., et al. 1999, MNRAS, 303, 659 
Houck, J.R., et al. 2004, ApJS, 154, 18 
Houck, J.R., et al. 2005, ApJ, 622, L105 
Hughes, D.H., et al. 1998, Nature, 394, 241 

Ivison, R.J., Smail, I., Le Borgne, J.-F., Blain, A.W., Kneib, J. -P., Bezecourt, J., Kerr, T.H., 
Davies, J.K. 1998, MNRAS, 298, 583 

Ivison, R.J., et al. 2002, MNRAS, 337, 1 

Ivison, R.J., et al. 2004, ApJS, 154, 124 

Jarvis, M. J., Rawlings, S., Eales, S., Blundell, K.M., Bunker, A. J., Croft, S., McLure, R.J., 
Willott, C.J. 2001, MNRAS, 326, 1585 

Kauffmann, C, Colberg, J.M., Diaferio, A., White, S.D.M. 1999, MNRAS, 303, 188 

Klaas, U., et al. 2001, A&A, 379, 823 

Knudsen, K.K., Kneib, J.P, Egami, E. 2006, |astro-ph/0604073| 



-22 - 



Kovacs, A., Chapman, S.C., Dowell, CD., Blain, A.W., Ivison, R.J., Smail, I., Phillips, T.G. 
2006, ApJ, 650, 592 

Krabbe, A., Boeker, T., Maiolino, R. 2001, ApJ, 557, 626 

Kreysa, E., et al. 1998, SPIE, 3357, 319 

Laurent, O., Mirabel, I.F., Charmandaris, V., Gallais, P., Madden, S.C., Sauvage, M., Vi- 
groux, L., Cesarsky, C. 2000, A&A, 359, 887 

Ledlow, M.J., Smail, L, Owen, F.N., Keel, W.C., Ivison, R.J., Morrison, G.E. 2002, ApJ, 
577, L79 

Lutz, D., Spoon, H.W.W, Rigopoulou, D., Moorwood, A.F.M., Genzel, R. 1998, ApJ, 505, 
L103 

Lutz, D., Maiolino, R., Spoon, H.W.W., Moorwood, A.F.M. 2004, A&A, 418, 465 

Lutz, D., Valiante, E., Sturm, E., Genzel, R., Tacconi, L.J., Lehnert, M.D., Sternberg, A., 
Baker, A.J. 2005, ApJ, 625, L83 

Lutz, D., Yan, L., Armus, L., Helou, G., Tacconi, L.J., Genzel, R., Baker, A.J. 2005, ApJ, 
632, L13 

Martmez-Sansigre, A., Rawlings, S., Lacy, M., Fadda, D., Marleau, F.R., Simpson, C, 
Willott, C.J., Jarvis, M.J. 2005, Nature, 436, 666 



Menendez-Delmestre, K., et al. 2007, ApJL, in press (astro-ph/0610915) 

Neri, R., et al. 2003, ApJ, 597, L113 

Persic, M., Salucci, P. 1992, MNRAS, 258, 14 

Pope, A., et al. 2006, MNRAS, 370, 1185 

Rigopoulou, D., Spoon, H.W.W., Genzel, R., Lutz, D., Moorwood, A.F.M., Tran, Q.D. 1999, 
AJ, 118, 2625 

Sanders, D.B., Soifer, B.T., Elias, J.H., Madore, B.F., Matthews, K., Neugebauer, G., Scov- 
ille, N.Z. 1988, ApJ, 325, 74 

Schweitzer, M., et al. 2006, ApJ, 649, 79 

Siebenmorgen, R., Haas, M., Kriigel, E., Schulz, B. 2005, A&A, 436, L5 



-23- 



Smail, L, Ivison, R.J., Blain, A.W. 1997, ApJ, 490, L5 

Smail, I., Ivison, R.J., Owen, F.N., Blain, A.W., Kneib, J.-P. 2000, ApJ, 528, 612 
Smail, I., Ivison, R.J., Blain, A.W., Kneib, J.-P. 2002, MNRAS, 331, 495 
Smail, L, Chapman, S.C., Blain, A.W., Ivison, R.J. 2004, AJ, 616, 71 



Smith, J.D.T, et al. 2006, ApJ, in press dastro-ph/ 0610913) 

Spoon, H.W.W., Keane, J.V., Tielens, A.G.G.M., Lutz,D., Moorwood, A.F.M. 2001, A&A, 
365, L353 

Spoon, H.W.W., et al. 2004, ApJS, 154, 184 

Sturm, E., Lutz, D., Tran, D., Feuchtgruber, H., Genzel, R., Kunze, D., Moorwood, A.F.M. , 
Thornley, M.D. 2000, A&A, 358, 481 

Sturm, E., et al. 2005, ApJ, 629, L21 

Swinbank, A.M., Smail, I., Chapman, S.C., Blain, A.W., Ivison, R.J., Keel, W.C. 2004, ApJ, 
617, 64 

Tacconi, L.J., et al. 2006, ApJ, 640, 228 
Tecza, M., et al. 2004, ApJ, 605, L109 
Teplitz, H.I., et al. 2005, AAS, 207, #52.08 

Thomas, D., Maraston, C, Bender, R., Mendes de Oliveria, C. 2005, ApJ, 621, 673 

Tran, Q.D., et al. 2001, ApJ, 552, 527 

Tremaine, S. et al. 2002, ApJ, 574, 740 

Veilleux, S., Kim, D.-C, Sanders, D.B. 1999, ApJ, 522, 113 

Voss, H., Bertoldi, F., Carilli, C, Owen, F., Lutz, D., Holdaway, M., Ledlow, M., Menten, 
K. 2006, A&A, 448, 823 

Weedman. D.W., Le Floc'h, E., Higdon, S.J.U., Higdon, J.L., Houck, J.R. 2006, ApJ, 638, 
613 

Yan, L., et al. 2005, ApJ, 628, 604 



This preprint was prepared with the AAS IATgX macros v5.2. 



Table 1. Summary of SMG observations 



Name RA DEC Error Ref. SL1 LL2 LL1 

J2000 J2000 " (7.4 - 14.5 /tin) (14.0 - 21.3 /an) (19.5 - 38.0 /tm) 

sxn. of cycles sxn. of cycles sxn. of cycles 



SMMJ00266+1708 (M12) 


00:26:34.10 


+17:08:33.7 





8 


F00 




120 


X 


15 


120 


X 


30 


SMMJ02399-0136 (L1/L2) 


02:39:51.87 


-01:35:58.8 





6 


G03 










120 


X 


30 


SMMJ09429+4659 (H8) 


09:42:53.42 


+46:59:54.5 






L02 




120 


X 


15 


120 


X 


30 


SMMJ09431+4700 (H7/H6) 


09:43:03.69 


+47:00:15.5 





3 


N03 










120 


X 


30 


SMMJ10519+5723 (LE 850.18) 


10:51:55.47 


+57:23:12.7 






102 










120 


X 


30 


SMMJ10521+5719 (LE 850.12) 


10:52:07.49 


+57:19:04.0 






102 










120 


X 


30 


MMJ120517-0743.1 


12:05:17.86 


-07:43:08.5 





4 


D04 




120 


X 


10 


120 


X 


30 


MMJ120539-0745.4 


12:05:39.47 


-07:45:27.0 





4 


D04 










120 


X 


30 


MMJ120546-0741.5 


12:05:46.59 


-07:41:34.3 





5 


D04 




120 


X 


11 


120 


X 


30 


MMJ154127+6615 


15:41:26.90 


+66:14:37.3 





1 


BOO 




120 


X 


15 


120 


X 


30 


MMJ154127+6616 


15:41:27.28 


+66:16:17.0 





1 


BOO 




120 


X 


15 


120 


X 


30 


SMMJ16369+4057 (N2 850.8) 


16:36:58.78 


+40:57:28.1 






102 


240 x 5 


120 


X 


10 








SMMJ16371+4053 (N2 1200.17) 


16:37:06.60 


+40:53:14.0 






G05 




120 


X 


15 


120 


X 


30 



I 

Note. — Names in brackets indicate aliases used in the literature. 



stated otherwise) 



Fraycr^t_alJ 



References. — In terferometric positions adopt ed for the IRS ob s ervati ons (from VLA 1.4 GHz if not 
iBertoldi etaU ll2000h . D04: PdBI mm po sition of iDannerbauer et al.l d2004h . F0 0: OVRO mm pos ition of 
PdBI mm po s ition of iGenzel et al] faooj) . G05: 1" from final CO position of iGreve et"aH l l2005h . 102: 

iLedlow et all tiOOi ). N03: PdBI mm position of iNeri et al. | 11200311 . For sources with no position error stated in the references 
assume < 0.5" from the interferometer configurations used. 



Ivison et al,l 



BOO: 
20001) . G03: 
(2002). L02: 
we 
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Table 2. Fit results. The best fit is in bold type (see the text for details). 



Name 




M82 




IRAS F00183 


-7111 


NGC 4418 


M82 + continuum 


X 2 


2 




X 2 


2 




x 2 


2 


Vz 


x 2 


2 


&z 


SMMJ00266+1708 


1.5 


2.73 


0.01 


3.1 


2.75 


0.02 


2.0 


2.66 


0.01 


1.4 


2.73 


0.02 


SMMJ02399-0136 


17 


2.80 


0.01 


12 


2.86 


0.01 


21 


2.72 


0.01 


1.1 


2.81 


0.02 


SMMJ09429+4659 


0.58 


2.38 


0.01 


1.5 


2.43 


0.03 


0.86 


2.33 


0.01 


0.59 


2.38 


0.02 


SMMJ09431+4700 


10 


3.33 


0.01 


3.2 


3.26 


0.04 


4.3 


3.20 


0.01 


0.99 


3.36 


0.02 


SMMJ10519+5723 


0.68 


2.67 


0.01 


1.2 


2.67 


0.01 


0.80 


2.64 


0.02 


0.75 


2.67 


0.02 


SMMJ10521+5719 


0.74 


2.69 


0.01 


1.3 


2.75 


0.02 


1.1 


2.73 


0.01 


0.63 


2.69 


0.02 


MMJ154127+6616 


0.55 


2.78 


0.01 


0.82 


2.79 


0.04 


0.60 


2.67 


0.02 


0.35 


2.79 


0.02 


SMMJ16369+4057 


0.60 


1.21 


0.01 


1.68 


1.17 


0.01 


0.66 


1.12 


0.01 


0.30 


1.21 


0.02 


SMMJ16371+4053 


1.1 


2.38 


0.01 


1.2 


2.48 


0.01 


1.3 


2.34 


0.01 


0.41 


2.38 


0.02 



Table 3. Properties of the SMG sample 



Name 


S850 (im 


Sl.4GHz 


K 


Opt. 


Magn. a z h 




SpAH7.7fjm d 


Scont.7.7 fJ.m S 


T d 


L m x 10 13 


Best fit 




mjy 




mag 


mag 






mjy 


mjy 


K 


L o 




olvllVlJUUZDO-rl lUo 


lO.Oltz.4 


C\A _|_1 cf 


zz.o 




1 


Z. 10 


U. ( OU 


U. loU 


OO. O 


n k 
U.o 




SMMJ02399-0136 


23.0±1.9 g 


526±10 h 


17.8 s 


R=21.2« 


2.5 2.80 


2.81 


1.12 


1.44 


53.4 


2.5 


M82+AGN 


SMMJ09429+4659 


4.9±1.5 i 


970±3.5J 


19.7J 


R=25.2J 


1.3 


2.38 


0.683 




83.0 


6.3 


M82 


SMMJ09431+4700 


io.5±i.8 J 


55±3.5 j 


20.2^ 


R=23.8 j 


1.2 3.35 


3.36 


0.607 


1.04 


41.6 


0.9 


M82+AGN 


SMMJ10519+5723 


4.5±1.3 k 


47±10 k 


>20.4 k 


I=24.6 k 


2.69 


2.67 


0.491 




42.0 


0.5 


M82 


SMMJ10521+5719 


6.2±1.6 k 


278±12 k 


>20.6 k 


1=22. 7 k 


2.69 


2.69 


0.370 


0.122 


61.6 


2.9 


M82+AGN 


MMJ120517-0743.1 


6.3±0.9' 


40±13 m 


22. 5 m 


R~25.4 m 




>3.6 












MMJ120539-0745.4 


6.3±1A 1 


55±13 m 


>22.7 m 


R>26.2 m 




>3.6 












MMJ120546-0741.5 


18.5±2.4 1 


42±13 m 


21. 9 m 


R>26.2 m 




>3.6 












MMJ154127+6615 


10.7±1.2' 


81±13 n 


>21.2P 


R>24.9P 




>3.6 












MMJ154127+6616 


14.6±1.8' 


67±13 n 


20.5 n 


R>24.5 n 




2.79 


0.472 


0.159 


35.0 


0.9 


M82 +AGN 


SMMJ16369+4057 


5.1±1.4 k 


74±29 k 


18. 2 k 


R=22.5 k 


1.19 


1.21 


0.493 


0.205 


27.5 


0.1 


M82+AGN 


SMMJ16371+4053 


11.2±2.9° 


74±23° 


19.21 


1=23.21 


2.38 


2.38 


0.519 


0.242 


34.3 


0.7 


M82+AGN 



a Adopted lensing magnification where applicable, magnification 1 is assumed otherwise. See text for references. The submm and radio fluxes and 

optical/near-infrared magnitudes listed here are not corrected for amplification, while the infrared luminosities include the magnification correction. 
Magnitudes are on the Vega system. 

b Redshift from previous measurements. See text for references. 

c Redshift from best fit to IRS mid-infrared spectrum (this work). Uncertainty 0.06 including effects of potential template mismatch (See Sect. 4 
for details). 

d Flux of the PAH 7.7/xm rest frame feature after continuum subtraction. 
e Flux of the continuum 7.7 /im rest frame. 

f Fraver et al. (2000) s Smail et al. (20021 h Smail et al. (20001 1 Cowie et al. (2002) 

J Jjedlow_et_aL (2002). Opt/NIR magnitudes for SMMJ09431+4700 are dominated by the H6 component while H7 is the dominant submm component 
(Tacconi et al. 20061 . 

k Ivison et al. (20021 1 Eales et al. (20031 m Dannerbauer et al. (20041 n Bertoldi et al. (2000) ° Chapman et al. (2005) 
P Dannerbauer (20041 i Smail et al. (20041 
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Fig. 1. — Spitzer IRS low resolution spectra (solid lines) of the sample galaxies. The detected 
sources are shown together with the best template fit (dotted and dashed line). Their 
redshifts are listed in Table [31 The last four spectra show the IRS non-detections. 
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5 




Fig. 2. — Histogram showing the redshif t distribution fo r our S MG sample (solid) and for 
the submillimeter-flux-limited sample of IChapman et al.l (120051 ) (cross-hatched), scaled for 
the maximum value. The dotted line at z — 2.3 indicates the median value obtained by 
Chapman et al.l (120051 ) and the dashed line at z = 2.79 indicates the median of our sample 



distribution. Our value is calculated assuming z > 3.6 for undetected sources. 
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Fig. 3. — Average IRS spectrum of all 9 detected SMGs, individually scaled to S222 = 
15mJy in the rest wavelength. Both the 6.2 /zm and the 7.7 //m PAH features are clearly 
visible. The spectrum is well fitted by the starburst-like spectrum of M82 plus a weak 
continuum. 
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Fig. 4. — Histogram showing the ratio of PAH 7.7 /mi peak flux density and rest frame 222 /im 
continuum flux density for our SMG sample (solid) and for eleven local ULIRGs (cross- 
hatched). In this measurement of the mid-to far-infrared SED, SMGs are very similar to the 
local ULIRG population. However, they show a lowe r value than the low lumi n osity starburst 
M82. The M82 point is based on the PAH d ata of iForster Schreiber et al.l (120031 ) and the 
far- infrared continuum of IColbert et al.l (119991 ). obtained in large and similar apertures. 



